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Under the constraints of HQET the equations of motion of heavy meson distribution amplitudes 
of definite geometric twist, using the knowledge of their off-cone structure, are reformulated as a 
set of algebraic equations. Together with equations due to various Dirac structures various rela- 
tions between the (sets of) independent two- and three particle distribution amplitudes of definite 
geometric twist are derived and presented using both the notion of (double) Mellin moments and 
re-summed non-local distribution amplitudes. Resolving these relations for the independent two- 

- - - , particle moments in terms of three-particle double moments we confirmed the representation of ^±\„ 

^ ! by Kawamura et al. (Phys. Lett. B 523 (2001) 111). 
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^ ■ I. INTRODUCTION 

^3 : 

The intention of the present work is to extend, within the framework of Heavy Quark Effective Theory 
' (HQET) [l|, 0, 0, 0], our understanding on heavy meson wave functions or, equivalently, Hght-cone (LC) 

distribution ampHtudes (DA). In a previous work [3| we introduced already two- and three particle DAs of 
i definiteoeomeiric twist for pseucoscalar heavy mesons and related them to the usual LCDAs (of dynamical 

Dh; twist) 0,0,1,3. In addition, we derived various relations connecting separately the two- and three-particle 

^ . LCDAs. Now we like to derive additional relations connecting two- and three-particle DAs mutually. Such 

Q ' relations, based on exact operator identities due to the quark equations of motion (EoM) and equalities 

f~i , between the Dirac structures of the DAs, have been considered first for light meson DAs, see e.g. Refs. (lol . 

[m, s-iid later on, taking into account the simplifications due to HQET constraints, for heavy meson DAs, 
I see e.g. Refs. [I,0,[l3,[l3,[l3|. These relations, as long as the DAs are not expHcitly known off the light-cone, 

^ ■ have to be solved as differential equations for the heavy meson LCDAs. But now, since from recent work 

' [3, [13, the decomposition of various non-local QCD tensor operators into operators of geometric twist 

lO . is expHcitly known also off-cone all necessary differentiations of the DAs can be performed directly and, 

■ therefore, we need not to solve any differential equation to get the wanted relations. 

, To be more specific, the two- and three-particle DAs are related to the vacuum-to-meson matrix elements 

lO ' of generic bi- and trilocal operators which, for the present consideration, are given by 

o : 

. {0\q{x)U{x,0)TK{0)\B{v)) and {0\q{x)U{x,i3x)F^^{-&x)x''U{^x,0)TKiO)\B{v)) , (1.1) 

o ■ 

^ I respectively, where hy{x) and q{x) are the heavy quark and light antiquark field, respectively, \B{v)) is 

•r-^ . the (pseudoscalar) i?-meson state of (fixed) momentum P = Mv, T = {l,jaA'^a[3,"y5,j5jaAl5'^af3} is a 

' generic Dirac matrix and U{kiX, K2x) = Vexp {—ig J^^ drx^ Afj{Tx)} is the usual path ordered phase factor 

■ ensuring manifest gauge invariance {g is the strong coupling parameter and 1 > > some parameter). 
■ " " ' Let us remember that in Fock-Schwinger gauge {xA = 0), where the phase factor equals unity, the gauge 

potential A^{x) is related to the field strength F^i, according to A^{x) = /g^di? -dx" Fa^i'dx). 

In accordance with the definition of usual meson LCDAs [l^ but, additionally, respecting HQET con- 
straints the S-meson LCDAs arise by parametrizing the off-cone matrix elements l|l.ip . e.g., 
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{0\q{x)rK{0)\B{v)) ^IC'^{v,x) du(pa(?/) e-'"(^*), (1.2) 

Jo 

'\q{x)F^,{dx)x''Th,{0)\B{v)) = /C?^(«,x) / VuTa{u^,U2) e-'("i+''"^)(^-), (1.3) 
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where J^u = /o^diti /Q"^dit2 and x = x — v{^{vx) — y/ (vx)"^ — x^v^) , = 0, defines some Ught-ray being 
related to a; by a fixed non-null subsidiary four-vector which may be identified with the S-meson's velocity. 
As indicated, the matrix elements (|1.2p and l|1.3p are represented by the Fourier transform of the LCDAs 

w.r.t. the variable xP and additionally parametrized by a kinematic factor ]C^^^'^{v,x) depending on the 
tensor structure {a} of the nonlocal LC operator. These factors are the basic kinematical structures (of scale 
dimension s w.r.t. xd) of the matrix elements which, in principle, can be read off from their parametrization 
w.r.t. leading LCDAs of dynamical twist since at leading order geometric and dynamical twist coincide by 
construction. Their explicit form has been introduced in Ref. [5|. The integration ranges result from the 
fact that, in the framework of non-local LC expansion [13, HH, these matrix elements are shown to be entire 
analytic functions in the variable xP [l^l whose support is restricted to [—1,1]. Additionally, due to the 
(anti) symmetry of the relevant QCD operators the integration range can be restricted to < it^ < l,i ~ 1,2. 

Conventionally, LCDAs are characterized by its dynamical twist which, roughly speaking, counts powers 
of M/Q for the various terms in the kinematic decomposition of the matrix elements of non-local QCD 
operators Here, instead we use the original, group theoretically founded definition of geometric twist, 
T = (scale) dimension d — (Lorentz) spin j, which has been introduced for local QCD operators in Ref. [l^l 
and generahzed to non-local tensor operators on the light-cone in Refs. [2l,[2B,[23|- The decomposition of the 
non-local tensor operators into operators of definite geometric twist leads to corresponding decompositions 
of the LCDAs 

With the aim of applying the EoM we need the twist decomposition of local as well as non-local QCD 
tensor operators also off the light-cone which has been studied in Refs. The decomposition of non- 

local operators O^^r} with given tensor structure {a} into an (infinite) sum of non-local tensor operators of 
definite twist r formally reads 

r r 

with the projection operators 'Pj^j^'^ d) explicitly given in Appendix [A] up to second order tensors (d is 
the off-cone generalization of the inner derivative on the light-cone) . 

Considering bilocal off-cone operators, the corresponding meson DAs of definite geometric twist t, gener- 
ically denoted by ip'"j\u), are introduced according to Ref. [l3| (cf. also Refs. (26l. [2E 



(0|O|.j(x,0)|i?(«)) =^7'Pj^'^'>(x,d)/cW?j(«,x) /'due-'"(^-)<^M(«). (1.5) 
In the explicit computations, instead of taking the DAs, their (double) MeUin moments will be taken 

^a\n = / dw U" ipaiu) , (1.6) 

Jo 

Tain,™ = /W<-"< T,(zii, 1^2) =^ T,|„(^)= ^™T,|„,„, (1.7) 

•^0 m=0 ^"^^ 



(r) 

respectively. In terms of Mellin moments ip^^;^^ the two-particle matrix elements decompose as follows, 

{0\O,^,{xMB{v)) = Y.f:KJ}tli-^^^ (1-8) 

r n— 

where '^{^^[^^^^(a;, d) are the corresponding local off-cone projection operators (AppendixE])- In case of three- 
particle matrix elements {0\O[cr}ix,'&x,0)\B{v)) the moments (/s^^j^ have to be replaced simply by T|^j^^^(i?). 
Luckily, they are not required off the light-cone. 

An essential property of the geometric twist decomposition is, that it is uniquely determined by corre- 
sponding irreducible tensor representations of the Lorentz group (25l . [26| . Hence this method works for any 
type of matrix elements of tensor operators creating a power series in x"^ only and is independent, whether 
the matrix elements are defined by time-ordered field products or otherwise. Any logarithmic dependence of 
the matrix elements on x^, either being determined by the light-cone expansion of propagators as in (3ll | or, 
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in terms of , through renormaHzation by using renormahzation group equations (2ll . [2^ will be implicitly 
contained in the wave functions resp. DAs ipa{u) and Ta(ui,U2) as already pointed out in Ref. [s^l- 

The paper is organized as follows. In Sect, [ill we review the operator relations resulting from the quark 
equations of motion, show the structure of the relevant matrix elements in the trace formalism and introduce 
some operator relations resulting from relations between Dirac structures. In Sect. IIIll we present the most 
general off-cone structure of the bilocal axial vector matrix element in terms of independent (geometric) twist 
DAs, determine the (single) derivation of the required two-particle off-cone DAs, determine the necessary 
three-particle DAs of definite geometric twist and relate them to the conventional ones, and write the 
independent structure elements in terms of geometric twist DAs. In Sect. IIVI all independent relations 
between two- and three-particle Mellin moments are derived. In Sect.[V]the corresponding nonlocal relations 
between the two- and three-particle DAs are determined. In Sec. |Vl] the connection with the work of 
Kawamura et al. j§| is carried out and their result concerning Mellin moments <i>±|„ of usual DAs is re- 
derived; some remarks concerning the transverse momentum dependence of these DAs are made. In Sec. I VIII 
some conclusions will be drawn which might be of phenomenological relevance for B-physics. In Appendix 
A we review the explicit twist decompositions of relevant (non)local tensor operators. 



II. EQUATIONS OF MOTION AND OPERATOR RELATIONS 

The relations we are interested in are obtained by applying the partial derivative 9^ to the bilocal operator 
q{x)TU{x,0)hy{0) and q{x)j^TU{x,0)hy{0), respectively, off the light-cone. Therefrom, the following well- 
known exact operator identities can be derived, cf. Ref. [30|, 

a^g(x)7^r/i,(o) = q{x)^rh4o) + i f dd d q{x)F^,{dx)x''-f^rhy{Q), (2.1) 

Jo 

v^'d^ g(x)r/i„(0) = q{x)TD^K{{)) + i / di9 (i? - 1) q{x)F^,{dx)x''v^TK{0) + v^'5l{q{x)TK{Q)}, (2.2) 

Jo 

where = — lA^ and 5j^{-] is the 'total' derivative defined by 

8l{q{x)TK{Q)} = ^q{x + y)^hy{y)\^^^. (2.3) 

On the one hand, because of the well-known constraints of HQET, these equations simplify and, on the 
other hand, there exist various relations between them for different F-structures. Namely, due to the EoM 
for the light (massless) and the heavy quark, 

q{x) $ = Q and {vD)K = , (2.4) 

the first terms on the RHS of Eqs. I|2.ip and l|2.2p vanish. Furthermore, using the heavy quark on-shell 
constraint, 

i)K = hy , (2.5) 
and the well-known identities for the 7-matrices, especially 

7 7 7 - [9 9 ~9 9 +9 9 j 7^ + le 757i^ , (2-6) 

7a 7/3 = 9ap ~ iCTa/3 , Cr„/3 = (i/2) [7^ , 7^3] , (2.7) 

the pseudoscalar (75) and skew tensor (ij^Cap) structures can be related to the axial vector structure (7570). 
Some of these relations have been already derived for the geometric twist LCD As in our previous work [5|| . 

(1) Representation of 2- and 3-particle DAs using the trace formalism 

Before proving in general that the knowledge of the axial vector structure is sufficient let us show this by 
using the definitions of the 2- and 3-particle LCDAs in the trace formalism (see e.g. 0, Q as well as 0). 

First of all we consider the three vacuum-to-meson matrix elements containing the bilocal quark-antiquark 
operator for F = 75, 757q, 75ia'a/3, usually being parametrized by the help of the DAs $± = ^±{vx,x^), 
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together with their relevant derivatives: 

{0\q{x)rK{0)\B{v)) 



Tr 



1 + 



d^{m{x)rTKmB{v))\, 
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(2.8) 
(2.9) 
(2.10) 

(2.11) 
(2.12) 
(2.13) 
(2.14) 



where 9^ = d/d{vx) + 2x^ d/dx^ and, therefore, we have {vd) <I> = d^/d{vx) + 2{vx) d^/dx^. 

Computing these traces it is easily seen that for the determination of Ai . . . A4 only the Dirac structure 
r = 757q! must be considered (arguments vx omitted): 



(vd) (0|g(x)757a/i.(0)|B(z;))|^^_ = ifeM [ia Ai - vUvi)A2] /{vx), 
^^.{Q\q{x)l''l^lc.K{Q)\B{v))\^^^ = ifeM A3 + f„(ra)A4] /(wi) , 



(2.15) 
(2.16) 



Analogously, the vacuum-to-meson matrix element containing a trilocal quark-antiquark-gluon operator 
is parametrized in terms of four three-particle LCDAs "i! A{vS:;'d), 4'y(wi;i?), Xa{vx;i!}) and Ya{vx;'&) with 
■d being restricted to < < 1 as follows (ol.ls^l: 



{0\q{x)F^,{{)i)iTh,{Q)\B{v)) = 

{0\q{x)v''F^,{§i)S;TK{0)\B{v)) = 
{0\q{x)-f^'F^,{§S;)S;TK{0)\B{v)) = 



IbM 
2 

IbM 



Tr<^ 75r 



i + i) 



{y^,$ ~ (w£)7m) 



1-. - 1-1 



{vx] "d) 



(wx; §) - x^,XA{vi] ■&) + -^YA{vi; i9) 



{vx) 



,1 + ^ 



Tr<{ 75r^[i? Qi{vS:;^) ~ {vi) e2{vx;^)] }, 
Tr<{ 75r^ ^^{vii; ^) + {vi) ^^{vi; -&)] 



(2.17) 
(2.18) 
(2.19) 



Again, computing these traces for the Dirac structure F = 757a (arguments vx and § omitted) one obtains: 

(0|Q(i)w^F^,(i^i)i^757a/i.(0)|S(w)) = fBM[i^Qi-v^{v3:)Q2] , (2.20) 



with 



||q(i)7''F^,(^?i)i^757a/i.(0)|B(t;)) = feM [i^Q:^ + v^{vx)Q4\ , 

61 = [^-.4 + Ya] {vi; d) , 
62= [i!A + XA]{vi;d), 
63= [i!A + 2^v + XA]{vi;d), 
ei = 2[i>A~^v] {vi]d). 



(2.21) 

(2.22) 
(2.23) 
(2.24) 
(2.25) 



Putting expressions l|2.15p and (|2.16p and (|2.2Qp and (|2.2ip into the equations of motion, l(2TT1) and ((2?2|) . 
one arrives at the four differential equations of Ref. Q connecting the 2-particle and 3-particle DAs. How- 
ever, this approach has the drawback that the ^-dependence of the distribution amplitudes, especially their 
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dependence on x^, is unknown. In Ref. this is partially circumvented by compensating d^+/dx^ through 
combining two of the equations, disregarding another one containing also d^-/dx^ and finally solving only 
two equations which are independent of these derivatives. Below, in Sect.[Vllwe show that, despite omitting 
part of the information, their result concerning the Mellin moments $±|„(a;^ = 0) was complete. There, 
we also discuss the attempts (isl . [l5| to solve the problem of transverse momentum dependence of the 
B-meson wave function, i.e., to take into account the x^-dependence of <i>-|-|„(a;^ ^ 0). 

(2) General proof of sufficiency of the axial vector structure 

It is our aim to circumvent these drawbacks by using the decomposition of the appropriately parametrized 
2- and 3-particle vacuum-to- meson matrix elements into off- cone DAs of well- defined geometric twist whose 
x-dependence is completely known. Introducing corresponding (double) Mellin moments the necessary differ- 
entiations can be easily performed. When restricting to the Hght-cone they simply lead to ordinary algebraic 
equations for the corresponding LCDAs of (even) twists r = 2, 4 and 6. 

In principle, it is possible to restrict these considerations on the axial vector structure. At first, Eq. (|2.15p 
is sufficient to determine Ai and A2. But, looking at Eq. I|2.16p one observes that, in order to compute A3 
and A4, it seems to be necessary to know the twist decomposition of the second stage tensor operator or, 
having in mind relation l|2.7p . of the pseudo scalar and skew tensor operator. However, as will be shown 
now, the corresponding matrix elements can be related to that of the axial vector operator. 

Let us now derive the corresponding operator relations in the case of bilocal operators. Using the identities 
((TgI) and ((27|) one gets: 

9(a;)7^757a 7/3/11,(0) gap q{x)^^^^^h^{0) - q{x)^^^-fzicrai3h^{Q) 

- - [5^5/ - g^gc" + gcpg^"] qixh^-fphM ~ iec/" g(x)7p/i,(0). (2.26) 

Truncating both parts with 9^ one obtains a relation for F = 7510-0/3: 

d,,q{x)j''-f5iaa(3hy{0) da q{x)j5-fphy{0) - dp q{x)jr,jahv{0) + isap'"^ dp q{x)jahv{0). (2.27) 

Multiplying both sides of this equation with v'^ leads to a relation for F ~ 7570: 

df. q{x)r757M0) = [vf^ga^ ~ v^g^^a - vag^^] g(x)757/3/i.(0) + ie/^'iV qi^hMQ) 

= (v9)g(a;)757a/i.(0)-9a 9(^)75 /i„ (0) + 9^ g(x)7^75/i„(0) + ie^^^^Wp 9^ g(x)7,/i,(0) . (2.28) 

However, multiplying subsequently with 1;" does not lead to another independent relation for F = 75. 
Concerning the operator expressions on the LHS of Eq. I|2.2p we may state also the following relation, 

q{x)'-iziaaphy{G) = q[x)-fz[va^p - vp"ia)K[0) - iea/3'"^ Vpq[x)-fchy{Q), (2.29) 

which is obtained by multiplying (|2.26p with v^. 

When taking vacuum-to-meson matrix elements of relations (|2.26p - l|2.29p the vector parts accompanying 
the e-tensor do not contribute since S-mesons are pseudo-scalar, {0\q{x)^ahv{Q)\B{v)) = 0; therefore we get 

dp (0|g(a;)7^757a/i..(0)|B(i;)} = {vd) {Q\q{x)-i5-fahM\B{v)) 

- (vadp + Vpda) m{x)^5l^K{Q)\B{v)) , (2.30) 

dp {0\q{x)YlzicjapK{Q)\B{v)) = da (0|g(.T)757/3/^«(0)|B(«)) - dp (0|g(a;)757a/i«(0)|-B(v)) , (2.31) 

{vd) {Q\q{x)-i^icTaph.M\B{v)) = {vd) {0\q{x)l5{valp " V p-ia)K{Q)\B{v)) , (2.32) 

as well as 

{vd) (0|g(x)75/i.(0)|B(iO) = va {vd) {{)\q{x)-i5l''h,{0)\B{v)) (2.33) 

for the scalar operator. Obviously, all these relations are valid only under the constraints of HQET. 

Prom these relations it becomes obvious that all the expressions Ai, . . . , A4 may be derived by the help 
of the first derivative of the axial vector DA alone, either using Eqs. (|2.15p and (|2.16p or, eventuelly, using 
the structure of the pseudo scalar and skew tensor matrix elements. Therefore, only the geometric twist 
decomposition of {Q\q{x)^^^phy{Q)\B{v)) has to be considered in detail. 

Concerning the tri- local operators and their matrix elements the situation is more difficult. At first, 
relations similar to l|2.26p - l|2.29p also hold for the tri-local operators q{x)Fpy{'dx)x^^'^Th^{0) and 
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q{x)Ffj^^{'dx)x'^Thy{0) on the RHS of Eqs. (|2.ip and l|2.2p . respectively. But now, taking matrix elements, 
the e-terms, when multiplied with the vector part, ie^ip'"^ {0\q{x)Fp,y{'dx)x'' ^ahv{0)\B{v)) , will contribute. 

Furthermore, looking at Eq. (|2.2ip we notice that only the twist decomposition of the axial vector LCDA 
is necessary to determine 63 and 64. However, the determination of Qi and 62 is much more involved. 
This becomes obvious considering the axial vector operator q{x)Ffj_^{'dx)x'^j5hu{0) which by using the trace 
formula (|2.17p differs in their structure from (|2.20p but is simply given by fsMxa [©i — 02]- Therefore, 
has to be considered as an external parameter being irrelevant for the twist decomposition which, for dimen- 
sional reasons, has to be parametrized by a;^ and not by v^l So, the computation of Oi and 02 requires the 
determination of the twist decomposition of the second stage tensor 3-particle LCDA which can be decom- 
posed in an antisymmetric and a symmetric part. Whereas the antisymmetric part is easily parametrized by 
V[aii3], the parametrization of the symmetric part is troublesome, because gaf}{vx),VaVp{vx),VaXi3 + XaVjs 
as well as XaX0/{vx) may occur. This will be undertaken in the next section. 



III. TWO- AND THREE-PARTICLE DISTRIBUTION AMPLITUDES OF DEFINITE 
GEOMETRIC TWIST AND RELEVANT DERIVATIVES ON THE LIGHT-CONE 



In this section we collect all ingredients of the geometric twist decompositions of relevant tensor operators 
which, due to Eqs. (|2.15p and (|2.16p and (|2.20p and (|2.2ip . are necessary as input in the equations of motion. 
This will be done in terms of Mellin moments. By the way, there are three different types of matrix elements 
we are concerned with, namely 

• bilocal matrix elements: we only need their off-cone twist decompositions up to terms proportional to x^, 
since - after carrying out a single derivative - all the higher order terms vanish when projected onto the 
light-cone (i^ = 0), 

• trilocal matrix elements: we only need their on-cone twist decomposition which for the structures F^^{-dx)x'^ 
as well as Ff^,y{-dx)^^x'^ is known partly from our previous work and finally, 

• total translation terms which require a special treatment. 

(1) Derivatives of the two-particle axial vector distribution amplitudes on the light-cone: 
According to the foregoing considerations, especially taking into account the relation (|2.30p . it is only neces- 
sary to consider the axial vector DA and restrict the derivatives to the light-cone. However, there exist two 
independent parametrizations by and by x^/{vx) which will be considered separately. 

Let us first present the the off-cone (geometric) twist decomposition of the axial vector distribution am- 
plitude, resulting from Eqs. l|A.7l lA.S"]) . in terms of Mellin moments according to the definition l|1.8p (with 
the convention (^) = if either A; > j or fc < 0) : 



r n— 



E 



(-iPx)" ^ 1 



-x^ V nl 



jl V4(ra)V {n - 2j)l 



X Va {n + I — j — k) 



,fe jn-j-ky. 
{n + 2-k)\ 



with 



Xq ^ — ^ ( — iPx^^ \ ^ 1 



vx '■ — ' n 

n — 1 



,2 \J 3 



E 



{-\f{n-\)\{n-\-3-k)\(3 



! j! \^{vxY ) {n-\- 2jy. (n - k)\ 

^Sf =0 for n<j.. 



f%n = a^rd =0 for n < j 



(n-2fc)^^^+f , (3.2) 

(3.3) 
(3.4) 



resulting from the restrictions of the j-summation and being compatible with the fc-summation up to .7 + 1. 

(2) 

Furthermore, it can be shown that, besides the terms Vaip\^^ which appear for j = at arbitrary n, the 
u-parametrized part of the axial vector DA depends only on the combinations ^^j^^*^^ — ^^ai^^ whereas. 
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besides of the terms </5^2|n ^'^^ J = 0) its x-parametrized part depends only on the combinations </5^2|^^'^'' 

'^A2\n'''- Furthermore, in both cases, there occur no contributions of odd twist r = 3 + 2j. 
The first derivative of the axial vector DA, when restricted to the light-cone, reads: 



dp {0\q{xh57ah.M\B{P)) 



oo 

+ E 

oo 



(-iPx)" / VgVp (2) 5a0{vx) + (71 - l){VaXi3 + VfjXg) 



(vx) 



(— iPa;)" XaX/s n — 2 
n\ {vx)^A{n+l) 



2{n+l) 



(vx)^ 



(2) (4) 
fAl\n ~ fAl\r. 



.(6) 
llln 



iPS)" (5Q,;3(t;2;) + (??, — l)xaV/3 (4) aiai/3 (— iPi)" n — 2 



(vx)^ 



(4) _ (6) 
^A2\n ^A2\n 



(3.5) 



Surprisingly, the w-parametrized part, being related to LCD A 'P^j'i^, is symmetric in (a/?), whereas the 

(t) 

x-parametrized part, being related to LCDA V?^2|n' contains also a non-symmetric term. 



Now, multiplying by w'^, from Eq. I|3.5p we obtain 

f 00 

(vd) {o\q{xh,jMomp))\ = i/sM 5] 



(-iPi)" 



2(i;i) n! V n -f 1 

^ ^ 71—1 ^ 



(n + 2)^(f)„ + n^^7„ 



(4) 



(-iPi)" / n- 1 



i [(n + 2) + n ^%_] - [n + _ 3) ^f^^ 



(6) 



(3.6) 



and, multiplying by — .g"'', we obtain, analogous to (|2.16p . 



- a^(0|g(.T)7'^75/i.(0)|P(P)) 



A3 A4 = 



1 (-iP5)^ 



^\ E 



(wi) ^ 

^ ^ n— 1 



{-¥'!,1|„ + (« + 3)^^%„}, (3.7) 



whereas, using relation l|2.3ip . we find (in terms of ^^J^^Wi o^lyO 



9^(0lg(a;)7^75ifT„^/i.(0)|P(P)) 



2 ^ 2^ ^ ^^^(n - 1) ^Zr.^ (3.8) 



(wi) 



(wa;)2 '■ — ' n 



where, here and in the following, we use the convention a[c(6/3] — ^{ciabp — apba). 

Let us now express A^, fc — 1, . . . , 4 through the moments </5^j[„, i = 1, 2 of the axial vector DA as they 
follow from Eqs. ((2A5l) and ^EM - ^M- 



Ai = - 



E 



(-ips)" r n - 1 



A(vx) n\ 1 n + 1 

^ ' n=2 



{n + 2)^f^.+n^%^ 



-[n^%,_^ + {n-2)^%^ 



"^A2|n + ("-2)(p^^^|„ 



^2 9^,,-^, E 

A3 



1 {-iPiy 



2(vx) ^ n\ 

^ ' n— 1 

1_ ^ (~iPi)" 



n+1 



{n + 2)^%^^+n^'2iln 



(4) 



(vx) ^ n\ 

^ ' n=2 



A4 = - 



1 ^ (-ipi)^ 



- V 

(vx) ^ 



n^^^), +2(n + l)^W 



(3.9) 
(3.10) 

(3.11) 

(3.12) 
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(2) The relevant three-particle light- cone distribution amplitudes : 

After having determined the relevant bilocal LCDAs we now proceed to look at the trilocal ones occuring in 
Eqs. (|2.20p and l|2.2ip . Since in that case we are deahng with matrix elements on the light-cone, we can make 
use of the results already obtained in but, in addition, we also take into account those parametrizations 
which have not been considered there. Unfortunately, the twist projection operators do not commute with 
the vector w^. Therefore, before multiplying with we have to determine the twist decomposition of the 
matrix elements {0\q{x)Ffj,^{'dx)x'''-f5Thy{0)\B{v)) . However, due to its structure the matrix element must 
vanish when multipHed with x^^. Due to this, we are confronted with only one type of DAs in the cases of 
pseudo scalar and axial vector matrix elements and three types of DAs in the case of second stage tensor 
matrix element so that we get: 



(0|g(£)i^^,(z?i)i^75/i.(0)|B(«))/(/BAf) 



{0\q{i.)F^,{i&£)£''jMO)\Biv))/{fBM) 



n=0 



[vx) 



(^) 



n=0 



,(5) 

- P\r. 



(3.13) 



' lt^a]\n+l 



'x^ t[:} 



oo 



n=0 

ir)ll^' 

a] |n+l 



]{-iPxy 

nl 



(3.14) 



- T2\i 



{vx) 



r{r) 
' T3|n 



) 



„(r) 
- Tl\i 



{-iPxY 



2 ^ 

n=0 



{v^^Xa - {vx)g^a) 



,(4) 



' Tlln 



^ T2\i 

„(4) 
- T2\n 



V,_L'Va'X- 

2{vx) 



„(r) 
- T3|v 



2 r 



„(r) 
- Tl\r 



{vx) 



"(4) 
- Tllri 



"(4) 
^ T2\n 



„(r) 
- T2\r 



(^) 



E 



{vx) 



„(4) 
- T2\n 



-T 



(6) 
T2\n 



(3.15) 



We should remark that, contrary to the convention used in here, in order to be able to compare 
with the results in we are forced to define the tri- local matrix elements without the overall factor i, 
cf. Eqs. (1.3.8) - (1.3.17). (As a convention formulas from Ref. Q are indicated by writing "I" in front of the 
number being cited.) Of course, this does not change any of the results obtained in since there we only 
considered relations between two- and between three-particle LCDAs separately. 

The twist decomposition of the pseudo scalar and axial vector case have been given to simplify later on some 
notations. In the axial vector case, instead of and T^j, , which have been introduced in as related 



to VuXn and 



respectively, we now introduced T 



Tlln 



T 



Al\n 



related to antisymmetric and symmetric kinematic coefficients v^^x^-^ and v 



Al\n A2\n 



'^T3\n related to Xf^Xa/{vx) because it seems to be more appropriate to define the DAs in accordance with 
the various structures of well-defined geometric twist. 

Let us remark that all terms in the first equalities of (|3.13p and (|3.15p which contain x'^ can be ignored 
because, also after twist decomposition, they have a common factor x^ and, hence, vanish on the light-cone. 

Concerning the structure of the second stage tensor l|3.15p we should point to the following: An antisym- 
metric second stage tensor contains, in principle, 4 different twist structures but three of them vanish 

identically for the single possible parametrization; on the light-cone the remaining one T^^^|^ is of twist-4. 
On the contrary, a symmetric second stage tensor contains, in principle, 9 different twist structures [l^ (for 
each of the various parametrizations!) whereby three, namely two of odd twist and one of even twist, vanish 



respectively, as well as 



identically in our case. Therefore, the various LCDAs T 



1, 2, 3 for the symmetric tensor case should 



8 



occur with an additional subindex. However, due to the necessary vanishing after multipUcation with a;^ 
and after restriction to the Hght-cone that distinction can be removed. As a result we obtain contributions 
of twist-4 and twist-6 which are the leading terms of two independent infinite towers whose additional higher 
twist terms appear with nontrivial powers of x"^ and hence vanish. 

In addition, we remark that in the above list of three-particle LCDAs the skew tensor structure is missing. 
The reason is that up to now we are unable to construct the twist projection operator onto a tensor of third 
stage having mixed symmetry, i.e. being neither totally symmetric nor totally antisymmetric (cf. Ref. @). 
However, due to the result of the preceding Section, we do not need it for the aim of the present paper. It 
would only act as a consistency check. 

Now, let us multiply with in order to get the twist decomposition of the relevant tri-local matrix 
elements: 



{0\q{x)v^F^,{di)i''j5hAQ)\B{v))/fBM - ivi){ei - 62) = (vi) LjI^t^^)j^) , 



(3.16) 



{0\q{x)vf'F^,i^i)i''-f5jM0)\B{v))/fBM - J2 



(-iP£-)Vl 



= [-^Tl|; 



+ T 



(4) 
T2\r. 



1 ~ 

n a;, 



i-iPiY' 



n + 1 



"f(4) 
T2\r 



-r(6) 



and, using the results of Ref. 



W , (3.17) 



||q(£)^^^,(i^i)7'^i^757a/i„(0)|S(«))//sA/ = v^ivi)J2 



1 ~ 

2 ^ 



n = l 



{-iPxY 



n+1 



(3) 



n 



(5) 



Al\n Al\n 



W ■ (3.18) 



Let us now express the 'dk,k = 1,...,4 through the moments Tti,* = 1,2,3 and flAiji — 1,2, of the 
axial vector DA as they follow from Eqs. I|2.20p and (|2.2ip [ observing T^g^, = = f^^2|o] ■ 



61 
62 
63 
64 



1 ^ i-iPxY 

2 ^ nl 

00 



„(4) 
- Tl\n 



„(4) 
- T2\n 



W + ( 2T^lW - 



n + 1 



^T2\n ^T2\n 



ji=0 



l^{jAPxY_ 



ri=0 
00 



71+1 



(3) 



(5) 



Al\n 



(19) 



n=0 



Comparing Eqs. ((3?T9|) - with Eqs. ([2:22)) - ^(TI^ we find 

t(4) 



^ _ 1 /-^(4) 
^ A\n 2 ' 



X 



A\r 



Tl\n ' ^T2\n 

(4) _ -^(4) 
Tl\n ^T2\n 

(6) n 



"^31" ' 2(n+l) [^T2\n - 

^ _ 1 /^(4) ^(4) _ (-,(3) 
^V\n — 2 \^-^Tl|n ^T2\n "Al|ri 

where the first three relations could be read off also by comparing Eq. I|3.15p with Eq. (I. A. 21) of Ref. 0; the 
last relation is obtained by observing 94|„ = f2^j|„ = 2'I'^|„ — 2^v\n- In we also showed '^v\n = '^vln 



-(4) 



"(6) 
- T2|i 



id) 



(3.19) 
(3.20) 
(3.21) 
(3.22) 



for 


n > 0, 


(3.23) 


for 


71 > 0, 


(3.24) 


for 


77, > 1 , 


(3.25) 


for 


71 > 0, 


(3.26) 
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and y4|„ - X^|„ — 'Ypjn 



T-iQi^ „ , ^ , I J- t^'tu- J- T~i--)U-. I " 2 I T'l|n. T^oi-^ I ~r J- r 



T (4) , t(4) _ o(3) _ 9^(4) 



|n + 


04|„ we find 




lor 




97) 
(i.ZI ) 


for 


n > 0, 


(3.28) 


for 


n > 1, 


(3.29) 



T3\n 2(n+l)V ^2|n T2|)iy 2V ^1|" T2|ny ^ P|n 

(4) 

Tl\n ' ^T2\n "Al|n ^ V\n 

i) „ o(5) ^ o(^) - 9T(4) 

Al\n 2(?7. + 1) V"^l|" Mlln^ ^ "A2|ri ~ ^-^Tl|n ' -^T2|i 

which allows, by using Ty|^ and Tp|^^, to circumvent the complicated expression for 6i|„ and 63|„, and 
also to avoid f2yii|„- 

By the foregoing considerations we observe that the four independent distribution amplitudes can be 
chosen as T^^i^, T~p|''„ ^^"^ '^v]n- This allows to write the (double) Mellin moments of &k as follows: 

eMnW = Ti;^lj^) + Tf^J^), (3.30) 

e^inW = r'^^l^i^) , (3.31) 

e3|„(^) = T^\VW + 2Tfv|nW' (3.32) 

QMnW = ^TUnW + - "^^vlW ■ (3.33) 

Let us mention that, according to the definition (|1.7p of the i?-moments, by comparing powers in i9, the 
relations 1(3:23)) - (|3^ hold equally well for the double moments T^^^ ^ itself. 

Finally, let us put together all those lower three-particle moments which vanish due to the foregoing 
considerations. According to the definitions l|3.13p - p.lSp as well as relations l|3.27p - p.29p we find 

for i = 2, 3 , (3.34) 
for i = 1, 2 , (3.35) 

(3.36) 
(3.37) 





^(6) 
■^Ti|0 


= 






= 


o'r(4) 




~ "ai|o' 


^(4) 
■^T1|0 


-^(4) 
T2|0 





(3) Remark on the total translation parts : 

Finally, there is the last term in (|2.2p . related to the total translation, to be considered. Yet, it is not known 
how to deal with these total derivatives by the method of twist decomposition. For that reason we introduce 
an "effective mass" A = M — mi, in style of the first of Refs. Q and simply define 

ivf^Sj; {{0\qix)TKmB{P))}l^^ = A {0\q{i.)TK{0)\B{P)) . (3.38) 

In the limit of infinite heavy quark masses A goes to zero. Furthermore, the contributions of the total deriva- 
tive v^Sj^ {{0\q{x)Thy{0)\B{P))}\^_^ may be assumed to be small. Furthermore, this substitution must be 
done only for the axial vector case, v'^S'-^{{0\q{x)j5jahv{0)\B{P))}\^_^ — > —lA {0\q{x)j5jahv{0)\B{P)) , 
whose twist decomposition already has been determined in the preceding Section. 



IV. RELATIONS CONNECTING TWO- AND THREE-PARTICLE MELLIN MOMENTS 

Now we are ready to discuss in detail the relations (|2.ip and i|2.2p for the (only relevant) axial vector 
case which connect the matrix elements l|2.16p with (|2.2ip and l|2.15p with (|2.20p in combination with the 
corresponding total translation term, respectively. Namely, from Eq. (|2.ip with F = 757^, by applying the 
EoM we obtain 

(0|9(x)7'^757a/i.(0)|B(«))|^^- = if (0|g(£)F^,(^?£)i"7'^757a/i.(0)|S(iO), (4.1) 

Jo 
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whereas from Eq. (|2.2p after applying the EoM and using relation l|3.38p we obtain 



(vd) (0|g(x)757c./i.(0)|B(«))|^^^ = i (0\q{i)F^A^i)i''v^j5laKmB{v)) 

Jo 

-iA (0|g(i)757a/i.(0)|B(P)). (4.2) 

Below we derive also a relation, connecting both the left hand sites of Eqs. (|4.ip and (|4.2p without using 
the total translation part. This leads to another relation connecting bi- and tri- local matrix elements. 

Due to the appearance of the derivative, the two-particle LCDAs contain an overall factor l/{vx) whereas 
the three-particle LCDAs contain an overall factor (vx) and the total translation terms due to l|3.38p do not 
have such an additional factor. Therefore, in order to derive the wanted relations, we have to compare equal 
powers of (vx) or, equivalently, powers of (— iPi). 

(1) To begin with, we consider relation (|4.ip by writing down the Va{vx)- and Xa- terms of relations (|2.16p 
and (|2.2ip separately: 



,^(-iPi) 



oo 



{-iPx) 



{-iPx) 



M ^ 



2(n + l) 



"All; 



(5) 
A2\i 



where the relative factor l/{vx)^ between left and right hand side has been rewritten as — M^/(— iPi)^. 
Comparing equal powers of (— iPi) and using relations (|3.28p and l|3.29p . we find (remember (p^||o = 0) 



^v\n-2)i^) forn>l, (4.3) 
-M^n 1) = n{n -1)1 ^ (t^^^^ + 2T^^)„_ J (^) forn > 2 . (4.4) 



The last equation can be used to replace <i5^2|n ^^^^ '^^^ Set 



dz9 19 ((3n + 1)T^^]|„_2 + (n - 1)T^"^|„_2 + 2(n + 3)t|;|^„„2) W for n > 2 . (4.5) 



(2) Next, we consider relation (|4.2p by writing down the Va{vx)- and Xa- terms of relations l|2.15p and (|2.20p 
separately: 



M2 ^ (-iPi)"-l 

— L — -I — 

n=l 



E 



n + 1 
(-iPi)"+i 



n=0 

M2 ^(-iPi)"-i fn-1 



M (i?-l)T(f)j^)-HAA/^ 



(4) 
2\n 



n=2 



n + 1 
(4) 



i^ + ^)'P%r.+^^Alln 



+ (" - 2) ~ 2 + (n - 2) 



HPiT (2) 

n! ^-^lln^ 



E 



(-iPi) 



d^ (^-1)1 T^^^V + T^T^m 



?1 + 1 



^(4) _ -^(6) 



(4.6) 



+ 24'L)W 



AA/^ 



(-iPi)" / 71 



n=l 



n! V 'T' + 1 



(2) _ (4) 



2^i1,„ 



(4.7) 
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Again, comparing equal powers of {~iPx) we obtain 



2{n+l) 



[n + 2) +n^%^ + '-^ ~^M^Al\n^i ^ [n - I) j M {i) - 1) T^^V-^W ^ (4- 



(2) 



{n - l){n + 2) - 2 n - (n + l)(n - 2) <^^°^| 



(4) 



2(n+ 1) 

- AM ( (?i - 1) 



(2) _ (4) 



= 2n{n — 



1 



+ - 2) ^^°^| J (4.9) 



(6) 



(n-1)/ d^(^-l)fT(f,V-2 + Tfp|U)W. 



,(5) 



where both relations hold for n >2, but the first one may be used also for n = 1. 

(3) In addition, we remark that from Eq. I|2.ip and relation (|2.28p together with the corresponding one for 
q{x)Ff^,^{i}x)x''^'^^5^ahy{0) = v'^q{x)F^,^{i}x)x''^'^"f5"fa7i3hviO) when using the identity ((2!6|) . we get 

{vd) {0\q{x)j5lahA0)\B{v))\,^^^ = v^do. {0\q{x)j,j^h,{Q)\B{v))\^^^ 

- [d^ (0|g(x)7'^75/i„(0)|B(«)) - dT? ^ (0|g(a;)F^,(,?x)x^7''75/i.(0)|S(z;)) 

But, due to Eq. I|2.ip with F = 75 the terms in the last line of that equation cancel each other and, using 
the on-shell constraint (|2.5p as well as l|4.ip . we obtain a relation connecting two- and three-particle DAs 
without the total translation part 

{vd) {0\q{xh5lcK{Q)\Biv))\^^^ = d,, {0\q{x)rj5jMQ)\B{v))\^^,^ + d^ (0|q-(x)75/i.(0)|B(z;))|^^. 

-if di9i9 {0\q{i)F^,{i}S:)i'' {g^^-f5-ie^^P^Vp-/^)K{0)\B{v)). (4.10) 
Jo 

The tri-local terms are given by Eqs. IpTTSl) and ifSTTII which introduce the DAs T^|',^(i9) and T^y^^Jd), 
respectively, the derivatives of the axial vector part (on the l.h.s.) and of the tensor part (first term on the 
r.h.s.) are given by l|2.16p and (|2.15p . respectively, whereas the derivative of the pseudoscalar part (second 
term on the r.h.s.) can be determined with the help of relation l|3.5p . eventually using the on-shell condition, 
as follows: 

a„ {o\qixh,h4o)\Biv))l^^ = i/^Af 1^ £ tlEp: 



n + 1 



in + 2)^%^ + n^'2i\n 



(4) 



2"^A2|, 



{~iPx)" fn - 1 



- [{n + 2)^%,^ + n^f^^^] - [n^%,^^ + {n-2)^%^ 



+ 2{n-2) 



(4) _ (6) 

^A2\n ^A2\r 



(4.11) 



Let us observe that, due to the symmetry in (a/J) for the terms V'Aiin in l|3-5p the last result differs from 

p.6p only in the terms V^aIiu- 

Now, putting together all these expressions we obtain, after cancellation of various terms and already 
omitting corresponding sums, the following very simple expressions: 

Xn 



= 
= 



M^{Tf^„-2T^^l){^), 
[n + (n + 3) J + 2 n(n - 1) ['m ^ T<^^1_^ (^) . 



(4.12) 
(4.13) 



Comparing the last equation with the sum of Eqs. I|4.3p and l|4.4p we find another relation between the 
trilocal DAs, namely 







^{2T%,^ + T^^l + 2T<^l)m. 



(4.14) 



12 



Obviously, this expression can be used to replace relation l|4.13p . Here we should remark that, contrary to 
"strong" relations f^l3^ - l(336l) and l(337)) - (029]), the relations l(4J2l) and (|414l) hold for the integrated 
DAs only. 

At this stage we should remark that no further relation can be derived. One might think that relation 
(|4.10p . together with (|4.2p . can be used to express the total derivation part as follows: 

v^{0\6l{q{x)j57M0)}\B{v))\^^^^ a^(0|g(x)7''757«/i.(0)|SM)|^^-+a„(0|g(x)75/i„(0)|i3(«))|^^_ 

-if {0\q{£)F^,{dS:)x-' (5^„75 - ie/''"«p7.) h,mB{v)) 

Jo 

-if d^{d-l) {Q\q{S:)F^,{di)i''v''-f5laK,{0)\B{v)) . (4.15) 
Jo 

However, due to relations l|4.12p and l|4.13p this expression falls back to the original EoM (|4.2p . And, on the 
other hand, when using Eq. I|2.28p together with l|4.ip according to 

v^{0\dj;{q{xh57M0)}\B{v))\^^^ = a„ (0|g(x)75/i„(0)|i3(z;))|^^_ 

+ if diSd (0|g(£)i^^,(i?i)i''7^75 (7« - v^) /i,(0)|S(«)) 
Jo 

-if Mid-1) {0\q{i)F^,idi)i''v^-f5-fM0)\Biv)) , (4.16) 
Jo 

having in mind relation l|3.38p . we are finally led to Eqs. I|4.8p and l|4.9p . 

The main results of this section are relations l(4!3|) - (|4?5|) . (|48|) and (|49| . (|4.12p and l|4.14p connecting the 
moments of two- to three-particle DAs of low twist. In addition, equations l|3.23p - l|3.26p . p.27p - p.29p 
and (|4.12p and (|4.14p connect the moments of various three-particle DAs. In the next section we convert 
these local results from the (double) Mellin moments into the corresponding ones of their nonlocal DAs. 

Let us emphasize that all these equations exhibit exact relations between the various contributions of 
geometric twist. Thereby, the leading contributions of twist r = 2 occur only for two-particle distribution 
amplitudes. Higher twist contributions (r = 4, 6) occur for both two- and three-particle distribution ampli- 
tudes. Thus, contributions of lowest twist are related to the lowest Fock state. These conclusions, of course, 
depend on the validity of the assumptions made by introducing the "effective mass" A. 

V. RELATIONS FOR TWO- AND THREE-PARTICLE DISTRIBUTION AMPLITUDES 

In this section we first reformulate the restrictions of two-particle DAs, especially the restrictions concern- 
ing their lower even moments, in terms of the non-local distribution amplitudes which can be interpreted as 
Burkhardt-Cottingham-like sum rules. Second, we reformulate the relations l|4.3p - l|4.5p . (|4.8p and (|4.9p . 
connecting two- to three-particle LCDAs for low twist in terms of the non-local light-cone distribution am- 
plitudes. And, finally, we reformulate the local relations l(333| - ^2^, (|3:27P - ((3?29l) and (|412l) and (|414| 
in terms of non-local three-particle LCDAs. 

(1) Burkhardt-Cottingham-like sum rules for distribution amplitudes of lower twists 

First of all, let us translate the restrictions which hold for the lowest Mellin moments of the two- and 
three-particle distribution amplitudes into their non-local form. 

From restrictions l|3.3p and l|3.4p for the two-particle DAs we immediately conclude the following: 

/ dwM" <^i\+^^'(u) = for < n < j , (5.1) 

Jo 

ip^^l{u)=0 and / dwM"(p^^2^^^^(u) = for 0<n<j. (5.2) 

Jo 

Analogously, for the three-particle DAs, according to Eqs. I|3.34p - p.37p . we obtain 

/ T^^^(wi,U2) = / I?uT^®^(ui,U2) =0, (5.3) 
Jo Jo 
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/ Vunf^{ui,U2) ^ / Vunf^{ui,U2)^0, (5.4) 
Jo Jo 

1 , , 

Vu {2T[^1 + Tf^ - nf^^ {uuU2) - , (5.5) 
)(ui,U2) = 0. (5.6) 



Vu (2Tf + T^^l - T^^^ 

Remind, that any odd twist DAs (p'"^^^'^\u) and T^j^^-'''(wi, U2) as well as any even twist DAs Q^^^^''\ui, U2) 
vanish. Obviously, these relations are Burkhardt-Cottingham-like sum rules for the axial vector distribution 
amplitudes. 

(2) Relations connecting two- and three-particle dispersion amplitudes of low twist 

The non-local version of the relations fTSl - (jUS]), fTSl and (|T9l) . I|4.12p and (|4.13p . connecting two- to 
three-particle LCDAs, cannot be reduced to the level of DAs itself but must be given in terms of integrations 
over the DAs. The reason is that these relations connect Mellin and double Mellin transforms. This is not 
the case for the lowest order of relations (|4.3p , (|4.8p and (|4.13p connecting only Mellin moments as well as 
(|3.27p . I|4.4p and l|4.12p connecting only double Mellin moments. 
Let us demonstrate the method of derivation on the generic case 



1 

— Va. = / fW ^n-sW for 



n > s . 



Multiplying both sides of that equation by (— iPa;)" /{n — s)!, using the integral representation 

dAA""'' for n>r, (5.7) 



1 



n — r + I 

together with the definition of Mellin (double) moments and summing up, we obtain 

g(-LPx)^ Z'duu" /dAA"-'-V(w)= f'ddf{{)) V ('^-P'^)" ' fvu{u^+du2T-'n{u,,U2). 

i^-^y- Jo Jo Jo t^s -^0 

Obviously, the summation on both sides results in exponential functions. The LHS can be rewritten by 
changing mA = A', exchanging A'- and u- integration and renaming the variables (A', u) — > (u, w). Analogously, 
the RHS can be rewritten by changing U2'& = 1?', exchanging 1?'- and U2-integration and renaming the 
variables {"d' ,U2) {u2,w). By this procedure we arrive at 

\u e-'^""- u'~- f'^w'-^iw)^ /'p7ie-n«i+"2)P. f' ^ f (!^) q^u,,w). (5.8) 

Ju W Jo Ju2 W \W J 

Let us mention that we need not restrict to the light-cone since the content of that formula does not depend 
on whether we write a; or i. 

In applying this to Eq. I|4.3p we divide it by n[n — 1), observe (n + l)l[n{n — 1)] = 2/(n — 1) — 1/n and 
finally get with s = 2, r = 2 resp. r = 1, and /(i?) = d: 

M'[\ue-'^P- f^wwU^liw) +2(2--) ^^^liw)] = /pue-'("i+"^)^-/'— ^f}i'](«i,H, (5.9) 

Jo Ju \ W/ ^ Jo Ju2 ^ ^ 

and for relation (|4.4p . taking into account Eq. p.29p . we obtain 

M^j\uue-'-P^ j\w ^f^{w) = j'^Vu ^ \^2Tf^ + T^^^ _ nf^ {u,,w), (5.10) 

where n'^^i = '^Ti + ^tI ~ 2T(^^ as given by Eq. (0211) • 
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Now, we consider relation l|4.8p . For n = 1 it relates only two-particle DAs and is almost trivial, 



M j^duu(^[?,^f,{u) + vf,{u) 



(5.11) 



whereas for n > 2 we find the following non-local relations, 



- u^) ip^'^^iw) + {up + u^) tp'-^^iw) + 2w{w - - AAf u; (p^^] (w) 



J a Ju2 w \w J 



(5.12) 



The non-local version of relation 114.91) reads 



du e^'^-f"^ 



,(4) 



2w{w~2u) ifi'^liw) ~ ifi'^^iw) 



MI 



Vai (w) - Vai (w) ~ 2w(p\l (w) 



^ Vue~ 

/ 



^ f dw U2 \ r (4), . , ^(5). X 

J no W \W y L ^ 



(5.13) 



(3) Relations between three-particle distribution amplitudes of low twist 

First, let us consider relations l|4.12p and (|4.14p . Taking into account the steps leading to RHS of Eq. (|5.8p . 
they may be non-locally rewritten as 



= / Vue 

p1 



If w 



{Tf~2Tf){u,,w), 



0= / Vue-'^^^+^-^P- /'^!^(2t(^) + t(^)+2TV"M(ui,«;) 

Ju2 W W ^ 



,(4) 



Second, concerning relation l|3.28p . the simplest of relations l|3.27p - p.29p . we get 



(3) 



Furthermore, the nonlocal version of l|3.27p may be obtained as follows: 



'^(4) _ ^(4) _ „-^(6) 



{ui,u2) + [i- I dxy 



^(4) _ ^(6) 
^ T2 ^T2 



1 .1 

dui / du2 e-''"i+''"='-f'^ 

"'0 

'dA f'du, /'d«2 e-'("i+*"^)^^ 
Jo Jo 



^ J- j5 'T" J- I ^2 ^ J- ^3 



(m1,M2) 
{U1,U2) , 

which, when integrated over "i? with an arbitrary function /(i?) as 



-f (4) „ r'-^) 

^ T2 T2 



Ju2 VW/ 



"r(5) I -Y-(4) 



-Y-(4) _ -y-(6) 9-»-(6) 
Tl -"^ T2 T3 



dw ^ /U2\ 



r(4) _ t(6) 

^ T2 T2 



(5.14) 
(5.15) 

(5.16) 



(5.17) 



{ui,w) 



(ui,w), (5.18) 
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may be used to replace |^Tp' + T^-^' j {ui,w) in relation l|5.13p . Furthermore, introducing new variables 

0<u = ui+U2<2 and A' = uX, then exchanging integrations over u and A' and renaming thereafter these 
variables as v and u, we may rewrite the r.h.s. of Eq. (|5.18p as follows, 



\ w / 



W \W/ V + U2 



dv 



-^(4) _ „^(6) 
^ T2 ^T2 



{v,w) , 



(5.19) 



where the support restriction < u < 1 has been used for the w-integration. Now, since is an arbitrary 
function and the exponentials build up a complete orthonormal system of functions we conclude that 







9"f (5) ^(4) -^(6) _ 9'v-(6) 





/■I dw 


-f(4) -»-(6) 
^ T2 ^ T2 









(V,M2) 



(5.20) 



Analogous to l|5.18p the z?-integrated relation (|3.29p can be reformulated non-locally as follows: 



w \ w 



(4) -^(4) _ ^(3) 



i TdA /'2?.e"'(-+-)^^- t—f(^ 

Jo Jo Ju2 ^ 



leading to the relation 







2T 



(4) ^(4) _ 1 /^O) 



17 



17 



(wi, W2) - ^ 



du 



17 



(5) 



A2l fiui,w) 

(ui,u-), (5.21) 
(w,M2). (5.22) 



Im V + U2 

For the special case of zeroth double moments, i.e. for n = 0, instead of l|5.19p and (|5.2ip one gets 

T-)„, "I 

P ' 2 



-^(5) 1 /-^(4) _ ^(4) 



,(6) 



T3 



(ui,U2) = 0, 



Pit 



2T^^i' + T^^^ - 17^^] - 17^'^] iui,u2) = , 



(5.23) 
(5.24) 



which, of course, is consistent with relations (|5.3p - (|5.6I 



VI. CONNECTION WITH THE WORK OF KAWAMURA ET AL. AND HUANG ET AL. 



This section is devoted to relate the two- and three-particle distribution amplitudes of definite geometric 
twist and their Mellin (double) momenta given above to the corresponding ones of Ref. and 

(1) Connection of the DAs ^±{vx, x'^) with the DAs V'^iin O'f^d ^^'22\n geometric twist r 

For our aim it is sufficient, to consider the axial vector part only and to do that up to terms of order 0{x'^). 

Taking into account the expressions l|3.ip and l|3.2p and definition (|2.8p . 



2{vx) 



we obtain $±(i'x,x^), as well as their difference, up to terms of order x^ as follows: 



[a>+-$_](«a;,a;2) = -^ 



{-\Pxy 



n=0 



n + l 



(2) (4) 
VAl\n - VAl\n 



^ ' 71 — 2 



71-1-1 



(2) _ (4) 
^Al\n ^Al\n 



(4) (6) 
fAl\n - Va1[ 



(4) (6) 
'^A2\n - '^A2\r. 



(6.1) 
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Jo 



(2) (4) 



w 



(2) (4) 

Vai ' Vai 



{w)-2^f^{u) 



dw {2w — 3w) ■ 



dw {w ~ 2u) 



(2) (4) 
Vai ~ VAX 



(4) (6) 

"Pax - Vai 



iw) 
M-2 



(4) (6) 
Va2 - 'PA2 



<^+{vx,x^) = - ^ 



(-iPx)" r (2) n(n - 1) 



n=0 
nl 



! I"^^!!" (ra)2 4(71+1) 



(2) (4) 



"PaiW + — ;— w 



(2) (4) 

Vai - Vai 



(-iPa:)" r 1 



n + 1 



(2) (4) 



A{vx)'' 



E 



(-iPx)" \2{n-l) 



(n+1) 
(n-2) 



(2) _ (4) 
'^Al\n VAl\n 



(4) (6) f, (4) (6) 



- I due-'^P"^ 







^ dw 

1/7 



(2) (4) 



/ dwe"'"^"^ / — w-2wW2m 



(2) (4) 

Vai - Vai 



(w) 



(4) {b , r, (4) (b) / N 



(6.2) 
(6.3) 
(6.4) 



(6.5) 



(6.6) 



where the sums go over the Melhn moments of $±|„(a;2) and the integrals over the distribution amphtudes 
$±(w;a;2). Thereby, it is necessary to have in mind that some of the lower moments of (pAi\n and (pA2\n 
vanish, related to corresponding Burkhardt-Cottingham sum rules. 

Let us first remark that the extra terms </3^2|n V^A2\n contribute only to <^-{vx, x"^) thereby indicating 
that they are higher twist contributions as will be natural for the subleading Furthermore, it is seen 
that by the additional (p^2|Ti^terms the light-cone DAs <^±(vx) are overdetermined. This was the reason to 
reject these terms in our previous work Concerning the x^-terms we observe that they are given in terms 
of differences — 'V^Ai\n^ ^'^'^ v'^A2\n ~ 'v'^A2\n^ which seems to be the case also for higher orders of x"^ . 

Again, these additional terms are overdetermined if the (pA2\n^ are taken into account, but would remain 
underdetermined and thus would lead to a difficult interplay between <i>+(ui) and $_(ui) if the V'A2|ri's are 
absent. In the following subsection we demonstrate that these terms necessarily occur. 

(2) Derivation of Kawamura's et al. result from our decomposition into DAs of definite geometric twist 
In Ref. [§| the x^-dependence of ^±\n has not been made explicit and implicitly taken into consideration 
only partly. Therefore, let us take into account only the ^^-independent part of $±|„ by solving the local 
relations (|4.4p , l|4.3p and l|4.8p - when the three-particle DAs are replaced by those of Ref. [§| - in that order 
with respect to the independent variables <p^2|n' ^'aiIu ^^'^ 'P^Ai\n ^'^^ n >2 a.s follows: 



/ d^ ^ [^A\n-2 + XA\n-2 + 24'^/|„_2] {^) 



5.7) 



M'v'-aL = 2 I d^ 



^{{n - 1) [*A|„_2 - *l/|„-2] W + {n+l) [*A|„-2 + Xa\„-2 + 2^v\n^2] w} , (6.8) 
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- (n + l)(n -1) I dd [^A\n-2 + ^A|„-2 + a^^^'yi^.-a] w}- (6-9) 



Relation l|4.9p could be resolved with respect to </5^j|„ + '^'^^A2\n using these expressions. But, taking the 
combination {n — 1) (|4.8p — n (|4.9p . we can avoid the use of expression l|4.8p . The result is (for n > 3): 

+ 2<^A2|n = ^^S^^^ " 1) t^^l-^ + ^.4|„-2] (^?) + 2n + y4|„-2] W} 

^ ^" " ^ ^ { [*A|„-2 - *y|„-2] (^) + 2(n - 2) A [vl;^|,^_3 _ vl/^|„_3] (^)} 

/ dT? 7? [Vl/^|„_2 + X^|„_2 + 2*v|„-2] (i^) • (6.10) 

Jo 



2 

M2(n - 2) 7o 



From relation l|6.5p we observe, that these twist-6 DAs occur in the a; -dependent part of $_|„ only, but 
in the combination <i5^j|„ — '^V^aIiu- Therefore, a further relation seems to be missing. Possibly they may 
appear if the full a;^ -dependence of all the two- and three-particle DAs is taken into account. 

For n = 2 the three-particle DAs are independent of ^ and simply given by the double Mellin moments 

"^Aio.o = ^'y|0,0 = 1>^H, Xa\q,o = 0, Ya\o,o = 0, (6.11) 

where X% and Xj^ are given by the chromoelectric and chromomagnetic fields in the B-meson rest frame [H]. 
Furthermore, according to (|3.4p for n = 0, as well as (|4.3p and (|4.8p for n = 1 the following holds: 

^aI\o = 0' + Mill = 0, U4^%, - 4-A^%, = 0. (6.12) 

Now, observing the normalization of ^±\o = 1 let us determine the lowest moments of ^'±|„ from the 
expressions l|6.3p and l|6.5p : 

(6.13) 
(6.14) 

i1|2 - fAVM^ + iAl,. (6.15) 

This is in full coincidence with the well-known result of as well as Eqs. (33) and (34) of Ref. 

Let us now consider the higher moments of ^±\n- First, $+|„ has to be determined by solving Eq. I|4.8p 
iteratively with the result 





-|0 = 


(2) 


= 1, 




-|o = 


(2) 


1, 


$H 


HI = 


(2) 

-Vai\i 


= |A/Af, 




-|i = 


1 (2) 
2'^-41|l 


= fA/M, 


$H 


H2 = 


(2) 
~VaI\2 


= 2AVAf2 + |a| H 




-|2 = 


1 (2) 
-3'^A1|2 


2 (4) o 
- 3'^A1|2 



n + 2lVA// ^ \M ^ \ £+1 

^ ' fe=0 ^ ^ £=0 ^ 

" '"^TTY *A|n-2-fc,f + (n + 1 - /c)A:A|„_2-fc,£ + (?1 + 2 - kfj^^V\n-2-kJ 



(6.16) 



where we used definition l|1.7p and performed the i?-integrals. The first term of that expression, not containing 
contributions from the three-particle DAs, by convention is called the Wandzura-Wilczek part (The 
notion "Wandzura-Wilczek part" is used quite differently in the Literature. Originally [33| it was introduced 
to denote that contribution to the twist-3 structure function 172 which appeared as "geometric" combination 
gY^i^) — ^9ii^)+ fx dy giiy)/y of the twist-2 structure function gi, like the combination (^^] (u)— (p^] (w) in 
Eq. (j6.2p appears when ignoring 2 ip^^l (u) ; for a more detailed discussion of distinguishing between "geometric" 
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and "dynamic" WW contributions, see e.g. Refs. |28 . l29l . l3a|). Remark also that there appears no term 

(A/M)" because the three-particle DAs start with n = 2. 
Next, we express $_|„ by plus an extra term as follows 



^-\r —'^4 



1 



£=0 



Up to a trivial change of indexes expressions l|6.16p and (|6.17p coincide with the final result of Kawamura et 
al. given by formulas (28) - (31) of Ref. ^. 
From this result we observe that 

• despite not exhausting all informations contained in the EoM the result of Ref. j§| is complete as far as 
the DAs $±|„ on the light-cone are concerned, 

• from the point of view of geometric twist $+|„ really is of minimal twist t = 2, but $_|„ contains also a 
(smaller) part of minimal twist r = 2 together with two different parts of higher twist t = 4, 

• the above derivation of Eqs. (|6.16p and l|6.17p shows the necessity of maintaining the contributions ip^^lln 
since otherwise these expressions would read quite different in contradistinction to Ref. [§| and, furthermore, 
due to Eq. we were led to J^M i? [^'a|«-2 + ^A|n-2 + 2'I'v|n-2] W = which together with (|i?T2l) and 
(|4.14p had strange consequences. 

(3) Some aspects concerning the -dependence of ^±{vx,x^) resp. transverse momentum dependence of 
their Mellin moments 

Looking at (|6.3p and l|6.5p we find that ^±(yx,x'^) for n > 2 contain a;^-dependent terms which are related 

to the differences <p^;^|„ — V^'a^u ^^"^ "^Ailn ~ "^Aiiii'* ^ 1,2 of DAs of consecutive twists. Of course, these 

x^-dependent terms for the Mellin moments lead to contributions of transverse momenta /cj^. 

In Refs. [l3 . omitting the contribution of 3-particle DAs, i.e. in "Wandzura-Wilczek- Approximation", 
the corresponding a;^-dependent equations for the Fourier transform with respect to the longitudinal sepa- 
ration t = vx (for convenience, we change to the variable w = Mw), 

$^^(a;,x2)=y' ^e'-*$^^(t,x2), (6.18) 

has been exactly solved whereby they found that both functions have a common x^-dependence: 

$^^(c^,x2) =0^^(c^)x(x2c^(2A-tj)) , (6.19) 

where (oj) = ^^["^t'^^ 9{lo)0{2A - w), (6.20) 

2A^ 

and x(a;^t^(2A-tj)) = Jo ( \x±\Juj{2A ~ oj) ) with = -x^ . (6.21) 



In principle, we would be in a position to confirm that result since we know from expressions p.ip and l|3.2p 
the full x^-dependence of x"^) but had to perform the Fourier transformation which is by no means 

simple. However, since the solution l|6.19p - (|6.2ip is an exact one we stop here. 

In Ref. [iBl the full problem has been tackled by assuming a common transverse momentum dependence 
also for each of the (Fourier transformed) 3-particle DAs. Despite being plausible that assumption has to be 
verified. In principle, the x^-dependence of the three-particle DAs can be determined in the same manner 
as we did it for the two-particle DAs by using the projections onto operators of definite geometric twist as 
given in the Appendix. This remains an open problem which must be postponed to another paper. 

In addition, some simplifying relations between the 3-particle DAs are introduced. And finally, a special 
model for the difference ^a{ui,U2) — 4'y(ui, U2) is required which, together with the other two requirements 
leads to the restriction Ya (1*1,1*2) = — ^'/i(tti, ^2) = Xa{ui,U2) and, from the view of the present work, 
seems to be very stringent. 



VII. SUMMARY AND CONCLUDING REMARKS 

With the aim of extending our previous work we used our knowledge about the explicit off-cone 
structure of QCD tensor operators of definite geometric twist (up to tensors of second stage) for rewriting 
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the (relevant) EoM - connecting the heavy mesons two- and three-particle DAs on the light-cone - into a 
set of algebraic equations for the (double) Mellin moments corresponding to these amplitudes. Thereby we 
have taken into account the heavy quark on-shell constraint and the well-known relations between Dirac's 
7-matrices in order to show that, in principle, it is sufficient to restrict to the axial vector structure T = 7570 . 

First, we found that two types of two-particle distribution amplitudes, ^p^J^l [for r = 2,4,6] and (p^J^ [for 
T = 4, 6], and five types of three-particle distribution amplitudes, T^-^\ [for r = 4, 6] and T^g as well as 
^Ai [for = 3, 5] and d^l 

occur; higher twists would appear if the EoM were considered off the light-cone. 
In comparison with our previous work we introduced two additional types of DAs. These sets of independent 
two- and three-particle DAs of definite twist are much larger than the commonly used sets of DAs consisting 
of $± as well as "if a, Xa and Ya, respectively. In Sect.nilby comparing the representation of the various 
matrix elements in terms of the corresponding DAs we derived the relations between conventional DAs and 
and those of definite geometric twist, Eqs. p.23p - l|3.26p . together with three relations, Eqs. (|3.27p - l|3.29p . 
connecting part of the DAs of definite geometric twist - some of them as special combinations - with two 
further sets, Ty^ and Tp^ . Thereby, we were able to select four appropriate combinations of three-particle 
DAs of definite twist, Eqs. p.30p - (|3.33p . which can be used to simplify the representation of the four EoM 
under study. In Sect. IVII the conventional two-particle DAs <i>± are represented (up to first order in x^) 
by the above mentioned DAs ^^Ai,i = 1,2, Eqs. (|6.3p - (|6.6p . showing that already on the Hght-cone the 
conventional ones contain three independent DAs of definite twist. 

Next, in Sect. IIV^ we presented four sets of algebraic equations, Eqs. I|4.3p and (|4.4p as well as (|4.8p and 

(|4.9p . connecting some combinations of the two-particle Mellin moments of definite twist, 'f^Ai\n^ ^ ~ ^' ^' '^ith 
the t?-integrated combinations p.30p - l|3.33p of (t9-dependent sums of) three-particle double Mellin moments 
of definite twist. In addition, two independent i?-integrated relations, Eqs. I|4.12p and l|4.14p . between the 
(double) Mellin moments are derived by combining the two kinds of EoM, Eqs. (|4.ip and (|4.2p . with the 
relation l|2.28p which was due to the Chisholm identity. In Sect. |V] these relations are reformulated non- 
locally as relations (|5.9p - l|5.15p between the two- and three-particle DAs. Furthermore, the relations p.27p 

- p.29p are reformulated non-locally in relations l|5.16p . I|5.20p and l|5.22p . Also in Sect. |V] we presented the 
vanishing of lower (double) Mellin momenta, Eqs. I|3.3p and l|3.4p as well as l|3.34p - p.37p . by corresponding 
Burkhardt-Cottingham-like sum rules Eqs. I|5.ip and (|5.2p as well as l|5.3p - (|5.6p . Obviously, the results of 
Sects. IV] and [V] are the main results of this paper. 

Finally, in Sect. IVH by resolving Eqs. (|4.3p . (|4.4p . I|4.8p and (|4.9p with respect to the independent two- 
particle MeUin moments ^^A2\n^ ^^Ai\n^ f^Ai\n v'Ai\n + ^"^^n' ^ (|6.10p . we were able, as a 
consistence check, to re-derive the result of Kawamura et al. [9| for the Mellin moments ^'±|„, given here 
by Eqs. I|6.16p and l|6.17p . In principle, the relations (|6.7p - (|6.10p can be used to express <i>±|„(x^) at least 
(partially) up to order by the three-particle double Mellin moments — if not a further relation were 
missing which would allow to separate 'f^Ai\n ^''A2\n- opinion this requires the consideration of 
the x'^-dependence in next order of all the participating two- and three-particle DAs of definite twist which, 
however, was not the aim of the present paper. For that reason we also were not able to further comment on 
the transverse momentum dependence of the two-particle DAs and their relation to the three-particle DAs. 

Despite being mathematically more complex, the group theoretically motivated use of the notion of ge- 
ometric twist allows for a very clear distinction between the contributions of different twist to the various 
matrix elements of physically relevant QCD operators and the corresponding DAs. Based on the quantum 
field theoretical framework, see Refs. (2ll . [2^ . it allows also for a different look at the conventionally in- 
troduced DAs, e.g., concerning the support of the various DAs or the appearance of new sum rules. Our 
study also showed that the twist structure of <I>_|„ is more complicated than usually assumed and requires an 
additional DA of twist-4; in addition, the usual three-particle DAs are shown to be appropriate combinations 
of three-particle DAs of definite geometric twist, thereby YA\n a nd A\n appear as complicated combinations 
of three different DAs each, compare expressions l|2.22p - (|2.25p with p.l9p - 113. 22p . cf. also relations p.23p 

- p.29p . These results bring some light onto the recently raised question [30| if the three-particle matrix 
element l|2.17p requires the introduction of additional DAs. 

Furthermore, from the point of view of the present paper, it seems to be possible to solve in that manner 
the problem of determining the transverse momentum dependence of the heavy mesons wave functions. 
Of course, concerning the two-particle DAs this requires the consequent use of Eqs. (|3.ip and l|3.2p and, 
additionally, to find the (infinite) off-cone decomposition of Eqs. p.l3p - (|3.15p and p.l6p - p.lSp : the latter 
becomes somewhat more difficult because of the appearance of additional independent twist structures. 
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APPENDIX A: OFF-CONE TENSOR OPERATORS OF DEFINITE GEOMETRIC TWIST 

Actually, we are interested in obtaining the twist decomposition of bilocal as well as of trilocal operators. 
As mentioned, only the tensorial structure of the operator is crucial for calculation and not whether it is 
created bilocally or trilocally. Hence, we deal with trilocal operators as if they were bilocal, but have to take 
care of the third field when writing the results. 

We denote the generic nonlocal off-cone (pseudo) scalar, (axial) vector and second rank tensor operators 
as follows: 

N{kiX,K2x), Oa{KiX,K2x), Map{niX , Hi2x) , 

where a possible pseudo structure is not labeled. The tensor operator of second rank splits up in an 
antisymmetric part M[Q^](Kia;, K2a:) and a symmetric part M(Q,^)(Kia;, K22;) from which also the trace 
M{kix, K2x) = Ma"{Kix, K2x) could be taken. The corresponding light-cone operators are obtained by 
replacing x ^ x and applying the constraint x^ ~ 0. 

In principle, we can perform the twist decomposition equivalently in the nonlocal representation as well 
as in the local representation. Here, we choose the local one. The relations between the nonlocal and the 
local operators are given by the Taylor expansion according to (restricting to ki = k; K2 = 0): 

7V(ra,0) = ^^7V„(x), Oa(KX,0) = J^O„|„(x), M^p{Kx,0) = '^M^p\^{x). (A.l) 

n— n— n— 

For the sake of a compact notation and also in order to make obvious the relation between the on-cone 
and off-cone version, we introduce the 'interior' differential operator which on the Hght-cone is given by [37| : 

d„/(£) = {(1 + xd)da - ^x^n}f{x)\ . (A.2) 

Its harmonic off-cone extension d and the complementary off-cone x-operator are given as follows: 

da = (1 + xd)da - ^XaO^ Xq ~ Xail + xd) - ^x'^da- (A-.3) 

The operators d^, x^, X = 1 + xd and Xap = xpda ~ Xgdg span the conformal algebra so(4, 2) as do the 
corresponding 'interior' operators on the light-cone, cf. Refs. [33|,[l3|- Especially, the off-cone operators obey 
also the following relations 

{X - l)dlaXf3] = - {X + l)xiadp], d(aX0) = X(^adl3) + Xgaf3 , (A. 4) 

Off the Hght-cone, the decomposition of a tensor operator of finite rank r with respect to the irreducible 
representations of the orthochronous Lorentz group is given by a finite series of traceless tensors having a 
well-defined symmetry type, cf., e.g. [s^. In our case, the possible symmetry types are restricted by the fact 
that any generic local operator Or\n{x) resulting from a Taylor expansion is a homogeneous polynomial of 
order n and completely symmetric w.r.t. the n indices which are truncated by the vectors x. However, this 
also allows the application of the polynomial technique [s^l- Consequently, these polynomials will vanish if 
more than n derivatives act on it. Therefore, the projectors exhibit an intrinsic termination which avoids 
the occurrence of undefined fractions or factorials (see below) [l7| . 

The twist decomposition of the relevant local tensor operators is given by Eqs. I|A.6P - (|A.10|) below. 
Thereby, the various contributions of twist r = (canonical) dimension — (Lorentz) spin are labeled by tq 
plus the higher order contribution due to the decomposition with respect to the irreducible representations of 
S0{3, 1). Thereby, tq is defined as the twist corresponding to the (fictitious) entirely symmetrized operator; 
in cases, where the operator can exhibit entire symmetry, tq also denotes the minimal twist of that operator 
and in cases, where entire symmetry is not allowed for that operator, as for the antisymmetric tensor of 
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second rank, we use tq as a counter only, cf . [I4I . The contributions resulting from subtractions of traces 
are numbered by j which count powers of (accompanied by powers of □) in the projector. 

All off-cone projection operators include the following projectors i?„ onto traceless homogeneous polyno- 
mials of degree n [3a |. 



k=0 



(-l)'=(n-fc)! , 2,fe^fc 



(A.5) 



which are sufficient for the formulation of the off-cone twist decomposition of scalar operators. For the vector 
and tensor case they contain, in addition, also specific tensor operators which are related to the symmetry 
type of the different twist contributions. The scalar projection operator has contributions related to even 
spin only whereas vector and skew-tensor projection operators exhibit even and odd contributions: 



{n+l-2j)\ 



Aij\{n+l-jy. 



(A.6) 
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{a0)\n V 4Jj!(n+l-j)! Ti I ; 
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(A.13) 



The finite series of scalar off-cone operators of definite twist l|A.6P have been given in Ref. [l^l using the 
harmonic extension of corresponding Hght-cone functions |37|. The off-cone vector operators of definite twist 
(|A.7p and (jA.sp have been determined for the first time in but are reformulated here appropriately; the 
expressions of even twist tq + 2j contain two different series of operators of definite twist with the second 
series starting at j = 1 . The off-cone antisymmetric tensor operators of definite twist l|A.9p and l|A.10p and 
the off-cone symmetric tensor operators of definite twist (|A.lip - (|A.13P in any dimension have been given 
for the first time in Ref. they are reformulated here appropriately for D = 2h = A dimensions. In the 
antisymmetric case the contributions of even and odd twist, tq + 1 + 2j and tq + 2 + 2j , contain two different 
series of operators of definite twist where, again, the second ones begin with j = I. In the symmetric case 
the contributions of even twist, tq + 2j and tq + 2 + 2j, contain three different and a single series of operators 
of definite twist, respectively, and the contributions of odd twist, tq + 1 + 2j, contain two different series of 
operators of definite twist. 

All the series of local operators of definite twist terminate at jmax — [■^^y^] with / being the number of 
free tensor indices T. Summing up all the twist contributions results in the non-decomposed local operator 
showing that the twist decomposition is a decomposition of unity into traceless tensors of definite symmetry 
type, i.e. irreducible representations of the Lorentz group. Obviously, the twist projection operators are 
applied to the non-decomposed operators: 

They obey the well-known properties of projection operators. 

Let us remark, that the reduction onto the light-cone obtains simply by using d — > d and a; — > i, = 0, 
resulting, especially, in Hn{x'^\0) = 1 and restriction to j = 0, 1. These on-cone operators have been given 
already in our previous work [3], Eqs. (I.E. 8) - (I.E. 19). Thereby, besides Eqs. (|A.4p . the following useful 
relations should be obeyed: 

x„i/„(x2|n) = i/„+i(a;2|n)a;„(l + a;9), {1 + xd)da, H^x'^lO) = Hn-i{x''\a) d^, 
x^dfi Hn{x^-\a) = Hn{x'^\a)xadf3, dc^xp H„{x:^\a) = H„{x^\a)do,xp. 

Here we presented the twist projectors thereby putting together different contributions of the same twist 
but of different symmetry type since we are interested in applying them for the twist decomposition of a 
distribution amplitudes which are related to their non-forward matrix elements in a given parametrization. 
From the group-theoretical point of view, such a summation losses some information about the structure of 
the distribution amplitudes which, however, could be re-covered, if necessary. For a more detailed discussion, 
see, Refs. 

The local operators 0'^l^{x) of definite twist can be summed up into non-local operators O'f \Kx,0) of 
definite twist by appropriately rewriting the various fractions in n which appear in Eqs. (KM - 
in terms of integrals thus arriving at the generic forms IjA.ip . This also works in the case of distribution 
amplitudes when determining it from their Mellin moments in Section Hill 
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